/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | foam-extend: Open Source CFD
   \\    /   O peration     | Version:     4.1
    \\  /    A nd           | Web:         http://www.foam-extend.org
     \\/     M anipulation  | For copyright notice see file Copyright
-------------------------------------------------------------------------------
License
	This file is part of foam-extend.

	foam-extend is free software: you can redistribute it and/or modify it
	under the terms of the GNU General Public License as published by the
	Free Software Foundation, either version 3 of the License, or (at your
	option) any later version.

	foam-extend is distributed in the hope that it will be useful, but
	WITHOUT ANY WARRANTY; without even the implied warranty of
	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
	General Public License for more details.

	You should have received a copy of the GNU General Public License
	along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.

\*---------------------------------------------------------------------------*/

#include "primitiveMesh.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

bool Foam::primitiveMesh::checkEdgeLength
(
	const bool report,
	const scalar reportLenSqr,
	labelHashSet* setPtr
) const
{
	const pointField& points = this->points();
	const faceList& faces = this->faces();

	scalar minLenSqr = sqr(GREAT);
	scalar maxLenSqr = -sqr(GREAT);

	labelHashSet smallEdgeSet(nPoints()/100);

	forAll(faces, faceI)
	{
		const face& f = faces[faceI];

		forAll(f, fp)
		{
			label fp1 = f.fcIndex(fp);

			scalar magSqrE = magSqr(points[f[fp]] - points[f[fp1]]);

			if (magSqrE < reportLenSqr)
			{
				smallEdgeSet.insert(f[fp]);
				smallEdgeSet.insert(f[fp1]);
			}

			minLenSqr = min(minLenSqr, magSqrE);
			maxLenSqr = max(maxLenSqr, magSqrE);
		}
	}

	reduce(minLenSqr, minOp<scalar>());
	reduce(maxLenSqr, maxOp<scalar>());

	label nSmall = smallEdgeSet.size();
	reduce(nSmall, sumOp<label>());

	if (setPtr)
	{
		setPtr->transfer(smallEdgeSet);
	}

	if (nSmall > 0)
	{
		if (report)
		{
			Info<< "   *Edges too small, min/max edge length = "
				<< sqrt(minLenSqr) << " " << sqrt(maxLenSqr)
				<< ", number too small: " << nSmall << endl;
		}

		return true;
	}
	else
	{
		if (report)
		{
			Info<< "    Min/max edge length = "
				<< sqrt(minLenSqr) << " " << sqrt(maxLenSqr)
				<< " OK." << endl;
		}

		return false;
	}
}


// ************************************************************************* //
